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I. INTRODUCTION 

The generation of coherent light at frequencies which 
are not easily reachable by lasers, and resonant but al- 
most lossless soliton-like light pulse propagation in a 
semiconductor are remarkable and useful optical phe- 
nomena. It has been recently showiP^ that the first ef- 
fect can be, in particular, realized by using the modula- 
tional instabilit}/^ (MI) of a long light pulse (continuous 
wave) propagating in a semiconductor and resulting in 
the parametric growth of equidistant spectral sidebands. 
The second effect, the so-called self-induced transmis- 
sion (SITm),^ has been demonstrated both numerically 
and experimentally for short intense light pulses. Both 
of these effects originate from the exchange Coulomb 
exciton-exciton interaction, which results in the nonlin- 
ear interplay of a semiconducting medium with a light 
pulse spectrally centered near the excitonic resonance. 

Like most of the nonlinear optical phenomena in semi- 
conductors, the mentioned effects can be theoretically in- 
vestigated by using the microscopic semiconductor Bloch 
equations (SBEs),^^ which however are rather compli- 
cated for both an analytical and numerical treatment, 
and require large computational facilities, being based on 
a many-body formalism. To obtain any relatively simple 
analytical results either the limiting behavior of SBEs^ ^ 
or a ppro ximate but macroscopic models of these equa- 
tiong^2Ii^ have to be considered. One of such macroscopic 
models, the so-called J-model, was proposed in Refs. El 
and [TT| where it was demonstrated that the model suf- 
ficiently well reproduced results obtained by numerical 
solution of the full set of microscopic SBEs for a certain 
range of incident light pulse intensities. Here we apply 
this model to study the mentioned optical effects. Unlike 
previous work,^ where we have investigated MI and soli- 
tary waves formation in semiconductors by using just a 
single nonlinear equation for the semiconductor polariza- 
tion, the present formalism describes also the dynamics 
of electron-hole density (or equivalently, of the so-called 



inversion) in a semiconductor and thus does not have the 
restrictions (i.e. weak excitation condition) of the model 
of Ref. 1, 

The structure of the present work is as follows. On 
the basis of a macroscopic analogue of SBEs developed 
within the J-model and coupled to Maxwell's equations 
for the electro-magnetic field of the incident pulse, we 
derive the system of governing equations (Sec.|llj). Then 
we use this system to analyze the stability of a long light 
pulse with regard to small perturbations. As a result 
we have found a strong MI of the pulse spectrum with 
a very large gain. All the results were also confirmed by 
a direct numerical simulation of the system of governing 
equations without additional approximations (Sec. |III[). 
To demonstrate the validity of the used model we have 
also applied it to describe the well studied SITm effect.^ 
After performing the slowly varying envelope approxima- 
tion (SVEA)^ both the analysis of the governing equa- 
tions and their numerical simulations revealed a good 
agreement with known theoreticaP^^ and experimen- 
tal^ results (Sec. IV). CdSe was selected as a representa- 



tive material in this work because it was a sem iconductor 
utilized in known experiments on SITm.^^^^ 



II. MACROSCOPIC SEMICONDUCTOR 
MAXWELL-BLOCH EQUATIONS 

It has been demonstrated in Refs. 9 and 11 that in 
the coherent regime, and for certain conditions, the com- 
plete system of microscopic can be sufficiently 
well approximated by a set of equations for macroscopic 
variables that are much more suitable for an analyti- 
cal treatment. Such macroscopic variables are the di- 
mensionless complex envelope p of the polarization field 
V{z^ t) = {p{z^ t) exp[z/c2; — iuot] + c.c.} /2 and the dimen- 
sionless inversion w{z^t) = 2A^eh(^,^) — 1 of a semicon- 
ductor, where A^eh is the normalized electron- hole density 
(Aeh ^ [0,1]).^ Here it is assumed that light is either 
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guided or polarized, and z is the longitudinal coordinate. 
In complex form this set of equations can be given as: 

dtV = -i[^P + Jp{w + 1) + nw], 

dtw = 2i[np'' -rf^p], ^ ^ 

where A = — uj defines the detuning of the 
pulse central frequency uj from the Is-exciton reso- 
nance cjx, ^ = dcvE/h is the Rabi frequency, while 
E is the complex envelope of electric field E{z^t) = 
{E{z, t) exp[ikz — iut] -\- c.c.} and dcv is the inter- 
band dipole matrix element. J is a phenomenological 
model parameter, which originates from the exchange 
(renormalization) terms of SBEs. In the low excitation 
regime, in which 1] <C cJb (^b is the Is-exciton binding 
frequency), J was obtained on a microscopic level and 
was defined in Ref. 9, Within the unscreened jellium 
model one finds J = 13a;b/3. We assume that the use 
of this definition of J in system ^ is well justified, at 
least, for the low excitation regime.^ This is the so-called 
J-model and it has been showri^ that it reproduces re- 
sults obtained by numerical solution of the full system 
of microscopic SBEs in the range of incident pulse inten- 
sities, which corresponds to the condition ^ < uj]j. A 
system similar to set ([T]) was also obtained from SBEs in 
Ref. [T2j However, here we prefer to use the formalism of 
the J-model, in which, unlike Ref. T2", the electric field 
is considered separately from the exchange (renormaliza- 
tion) terms of SBEs, which in their turn are included in 
the model parameter.^ 

System ([T]) is more general than the nonlinear equa- 
tion for the macroscopic polarization field previously used 
in Ref. []j to investigate nonlinear optical properties of 
semiconductors. Set ([T|) can be reduced to the previous 
case without inversion""" in the low excitation regime un- 
der the well-known approximation w —1 + 2|pp. For 
J = set ([T]) is equivalent to the optical Bloch equations 
(OBEs),^^ which govern nonlinear light propagation in 
systems of independent two- level atoms. A more detailed 
discussion on a limiting behavior of system ([T]) solutions 
will be given in Sec. |iy[ 

Obviously system is coupled to the Maxwell equa- 
tions for the electric and magnetic fields, 1-L{z^t) = 
{H{z^ t) exp[z/cz — iuot] + c.c.} /2. Written for envelopes 
without any additional approximations and in a proper 



dimensionless form together with set ([T]) they are the 
governing equations in this work: 

{dx^ = {'i^' — St)('0 + Ap) — ik'r]^ 
drp = - nL)P + J'P{^ + 1) + 

drw = 2i[iljp* - ijj^p] - 7;(^ + 1). 

In the literature such systems are referred as the semi- 
conductor Maxwell-Bloch equations (SMBEs). Here we 
have performed the next scalings and redefinitions: rj = 
H/Hq^ ip = Qto^ X = z/zq^ T = t/to^ uo' = cjto, k' = /czq, 
A' = Ato, J' = Jto, zo = cto/n, Hq = hn/{dcvto), 
X = dto, 7^ = 7x^0 7 7^ = Tr^o, c is the velocity of light in 
a vacuum, n is the non-resonant background refractive 
index, d = 2dl^/{7raQheoeB) is the photon-exciton cou- 
pling parameter, which naturally appears in Maxwell's 
equations after transition to the indicated dimensionless 
variables, £b = is the bulk background dielectric con- 
stant, £o is the vacuum permittivity, and ao is the Is- 
exciton Bohr radius. Parameter d defines the width of 
forbidden frequencies region in the dispersion relation 
of a photon-exciton coupled state.^ The exciton damp- 
ing 7x and relaxation 7r parameters are introduced here 
phenomenologically, and to is an arbitrary time scaling 
parameter, which will be properly chosen in the following 
sections. 



III. MODULATIONAL INSTABILITIES 

In analogy with Ref. IT) using SMBEs ^ we per- 
form the MI analysiJ^ - we analyze the linear stabil- 
ity of a long light pulse propagating in a semicon- 
ductor with respect to small perturbations. For this 
purpose we substitute the perturbed fields and po- 
larization envelopes 7^;p}(x, T) = ({'^o;^o;Po} + 
{a; g;b}{x,T)) ex.p[iqx] as well as the perturbed inver- 
sion w{x^T) = Wo -\- d{x^T) into system After that 
we define the small perturbations as {a;^;6}(x,T) = 
{ai]gi] hi} exp[mx — i5T] + {a2] 92] ^2} exp[zJT — 
while d{x,T) = 2doRe exp[mx — iST] because the inver- 
sion is a real function. Then, requiring that system ([2| 
is solvable, we obtain within the first-order perturbation 
theor}/^ the dispersion relation hz{S) for perturbations: 
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where q = uj^[l — Xwo/{^' + J'{wo + 1))]^/^^ ig the per- frequency. By using the relation between the steady-state 
turbation wavenumber, and 5 is the relative perturbation amplitudes po = — V^o'^o/(A' + J\wo + 1)) and a conser- 
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FIG. 1. The MI gain G{S) spectrum according to Eq. (jsj for 
^0 = 5, A' = 1.5, and 7^ = 0.267. (a) For J' = 65.4 and 
Wo > —1. (b) The same wq as in (a) but smaller J' = 6.54. 
(c) For J' = 16.35 and wq > —1. (d) The same J' as in (c) 
but different wo < 1. Insets at (a) and (d) show the same 
gain spectra but for a smaller vertical scale, u is positioned 
at the origin. 



vat ion law K;^+4|pp = 1, which follows from ([T]), one 
can obtain a 4-th order algebraic equation for wq'. 

{wl - 1)1 A' + J\wo + 1)1' + 4^olV^o|' = 0. (4) 

This defines wq for given ipo and system parameters. 
Among its solutions we select the real ones (either two 
or four): wq G (— 1;1). In the most of the cases there 
are only two real solutions - a negative and a positive 
one, and they correspond to a semiconductor, which was 
relatively weakly or strongly excited before the incident 
pulse arrival. However, of course, for certain combina- 
tions of parameters there exist four different real solu- 
tions, each of which corresponds to a different excitation 
regime. Therefore, setting certain incident pulse inten- 
sity and detuning, and choosing one of the mentioned 
excitation regimes, we completely determine the disper- 
sion relation (|3| as all the other parameters are defined 
by a medium. It is important to note, that, like in Ref. 1, 
the excitonic damping is naturally taken into account in 
Eq. (|3| by a proper final redefinition A' A'—ij^. After 
this the steady-state amplitudes and parameter q also be- 
come complex, resulting in the appearance of conjugated 
terms in Eq. ^ after the MI analysis. 

In the frequency regions, for which Im K:{6) < 
is valid, the exponential growth of perturbations 
takes place, defining the MI gain spectrum G{S) = 
2 max |Im hi{S)\. Since Eq. (|3| is of 4-th order in among 
its solutions we select, for a given 5, the one with the max- 
imum absolute value of the gain. The MI gain spectra 
for different J' are given at Fig. [T] where high resonant 
gain peaks clearly manifest themselves. The nature of 
these resonant peaks, which are strongly affected by the 
precise value of 7^, is discussed in Ref. [Ij A possibility 
to vary J' due to screening effects is shown in Sec. |IV[ 
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FIG. 2. (a) Numerically simulated long pulse spectral evolu- 
tion according to system ([2]), the initial pulse is = 0, T) = 
V^oexp [-(2T/Tw)^'^], ipo = 3, Tw = 200, m = 20, while the 
other parameters are: A^ = 1.5, J' = 1.96, and 7^ = 0.05. 
(b) The corresponding analytical prediction on the MI gain 
maxima positions by Eq. Inset at (b) shows the same 

gain spectrum but for a smaller scale of G. uj' is positioned 
at the origin. 



while here we only demonstrate how the change of J' 
can influence the MI gain spectrum (compare Figs. [TJa) 
andjljb)). Increase of J' can lead to a shift of MI gain 
peaks farther from the incident pulse central frequency 
uj' and to their slight gain enhancement (Fig. [TJa)). At 
Figs. [TJc) and [TJd) for a fixed J' (and other parame- 
ters) we compare the MI gain spectra in different excita- 
tion regimes. In the regime of pre-excited semiconductor 
{wq < 1, Fig. [TJd)) the MI gain peaks are shifted essen- 
tially farther from cu^ than those in the case of initially 
unexcited medium {wq ^ —1, Fig. [ijc)). However, de- 
pending on the selected parameters the opposite behavior 
is also possible. It is also worth mentioning that in the 
case when all the four solutions of Eq. Q are real, one 
of them usually corresponds to a regime with zero gain - 
the pulse is stable with respect to small perturbations. 

To confirm results of the mentioned analysis we have 
also performed direct numerical simulations of the full 
system ^ by using a 4-th order Runge-Kutta algorithm 
with the following initial conditions: iIj{x = 0,T) is a 
wide super- Gaussian pulse, r]{x = 0,T) = iIj{x = 0,T), 
p{x = 0,T = 0) = 0, and w{x = 0,T = 0) = -1, what 
means that the medium is not excited before the field 
arrival (see Fig.[2|a)). A comparison with the respective 
analytical prediction on the MI gain maxima positions in- 
dicates an excellent agreement (Fig.|2|b)). The regime of 
a preliminarily excited medium was also realized numer- 
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ically by using the initial condition w{x = O^T = 0) < 1 
and demonstrated qualitatively similar spectral evolution 
in accordance with the analytical prediction on the side- 
bands positions. At this point for calculations we should 
specify the representative material, which in this work is 
CdSe. For this semiconductor at low temperatures (e.g., 
8 K in Ref. 5, 2 K in Ref. 14), at which excitonic ef- 
fects are essential:^ hw^ = 1835 meV, fkJb = 15 meV, 
hd = 1 meV, fi7x = 0.265 meV, 7r < 7x according to 
Ref. fT4| ao o:^ 5.32 nm, dcv/e ^ 0.343 nm, and eb = 9. 
Now we choose to = 1/a 0.67 ps and for simulations re- 
move a large number uj' uj^/d from system ([2| by mul- 
tiplying its first two equations by \/uj' and performing a 
space rescaling zq zqX/u^ ~ c/{nuj) 0.04 /im. In ac- 
cordance with this scaling one can estimate, like in Ref. 1, 
from Fig. [I] the semiconductor MI gain ~ 100 /im~^, 
which is huge in comparison with that of, for example, 
silica optical fibers (~ 10 km~^).^ 

It is interesting to note, that a MI gain spectrum qual- 
itatively similar to the considered ones can be obtained 
even for J = (the case of OBEs),!^ something that was 
impossible in the model of Ref. [1] 



IV. SELF-INDUCED TRANSMISSION 
SOLITONS 

To demonstrate the validity of the used macroscopic 
model we now apply it to a well studied nonlinear op- 
tical effect. The SITm effect in semiconductors^ is an 
analogue of the known self-induced transparency (SIT) 
effect in idealized systems of two- level atoms .I^^IMi 
though both of them consist in almost lossless soliton-like 
propagation of short (subpicosecond in semiconductors) 
light pulses in a resonant medium, they essentially differ, 
what is discussed in details, for example, in Refs. [SHSl 
Here we only briefiy remind that due to strong infiuence 
of Coulomb exchange interactions between the Wannier 
excitons in semiconductors, which are not inherent to 
electrons in atomic systems, the Rabi frequency of the 
incident field approximately doubles^^ i n a se miconduct- 
ing medium, when it is about Q ^ ^^| iQ | ii | in particu- 
lar, this results in a substantial deviation from the so- 
called area theorem, which is known from the theory of 
a light pulse propagation in systems of two-level atoms 
and mainly states that only for incident pulse of an area 
= 2 \ilj\dT\x^Q equal to 27r SIT effect takes place 
- pulse is stable during the propagation. A pulse with 
an area, which is integer multiple of 27r, should undergo 
a breakup into 27r-pulses.^^ Such phenomena have been 
also found in semiconductors but fo r puls es with area 
TT and their respective multiples. ^^^^^ Such experi- 
mental and numerical observations (based on microscopic 
SBEs) have never been explained in clear physical terms 
in the literature - from which one of the motivations of 
the present work originates. 

Here we demonstrate the latter effect using system ^ , 
changing the variable r = T — x/V {V is the dimen- 



sionless pulse group velocity), and applying the SVEA.I^ 
SMBEs ([2| transform in this way into system: 



drP=-i[{A'- 
drW = 2z[V^J9* 



-V^*p]-7;(^ + i), 



(5) 



where we set to = ^/a^ and a = 2{duj)~^{l/V — 1) > as 
V « 1.^ Following Ref. UHl where light propagation in 
a dense medium of interacting dipoles has been investi- 
gated by using a system analogous to set ([5|, we briefiy 
reproduce an exact solution of Eqs. ([5| in the coherent 
limit (7x,7r = 0). Evidently, when using the first equa- 
tion in set ([5| the last one becomes integrable, and hence 
w = —1 -\- 2p^p. Therefore, the whole set ([5| reduces to 
a single equation for ip only: 

d^^ + i{A' + 2f\^\^)dr^ + {2\^f - 1)^ = 0. 

Looking for its solution in the form ^p{r) = 
?/^(r) exp[z(/p(r)], ip = |?/;|, we finally come to a couple 
of equations: 



d^ijj + Aip^ + Bijj^ - CV^ = 0, 



AV 



~2 



(6a) 



(6b) 



where A = 3J'74, B = AV' + 2, and C = 1 - A'V4. 



Eq. (6a) contains the so-called cubic-quintic nonlinearity, 
it is integrable and has a solitonic solution!!^ 



2Vc7b 



(l + v/1 + 16AC/(3B2) coshpyCr]) 



1/2- 



(7) 



The parameter space for this soliton is |A'| < 2, A' > 
—21 J' - it can form only if the incident pulse is spec- 
trally centered very close to the excitonic resonance. The 
area of soliton ([7|) 6>j = 2 dr tends to 27r only if 

J' ^ and VCsech[\^r]^ otherwise almost for 

the whole range of parameters Oj < 2i\ (see Fig. [3|. 
This, in accordance with known theories^ and exper- 
iments,^ ^ clearly demonstrates within the used model 
that exchange Coulomb interactions between excitons, 
from which parameter J actually originates,^ can re- 
duce the area of a SITm soliton in a semiconductor (see 
Fig. [sj^d), an explanation on why J can vary is given 
further in this section). 

To support the above mentioned analytical conclu- 
sions we have also performed direct numerical simula- 
tions of system ([5| after returning to the x, T- variables, 
and by employing a n umerical algorithm similar to that 
described in Sec. |III[ In Fig. Qa) a propagating pulse 
of a doubled area of dj o:^ 2.27r [parameters are given 
at Fig. [sj^a)] undergoes a breakup in two well separated 
pulses of areas Oj 1.27r and :^ Itt [see Fig.jljb)], what is 




FIG. 3. (a) Temporal shapes of IV^I^ (soUd), I^J^ (dashed), and 
w (dot-dashed hne) in accordance with Eq. ( TI for = 1 and 
= 3. (b) The same functions but for = and J' — 2. 
(c) The same but for A^ = — 1 and J' — 1. (d) Dependence 
of the pulse area Oj on the value of for A^ = { — 1; 0; 1} 
- dashed, solid, and dot-dashed lines, respectively. Grey dots 
indicate the pulse areas for the considered cases. The param- 
eter space restrictions of soliton ([7| are clearly demonstrated 
by the example of A^ = — 1. 



in good accordance with the prediction on a SITm soliton 
shape and area I.Itt given by Eq. ([7| (compare with 
Figs.[3]^a) andjsjd)). We should also notice that the soli- 
ton shown at Fig. IsFc) (A' = —1, J' = 1) was unstable 
under the propagation. 

Now, in order to perform a more detailed limiting 
analysis of solution ^ and to obtain some quantita- 
tive estimates we return to the microscopic definition 
of model parameter J, which, as we have already men- 
tioned, is given in Ref. ^ ( J = /3i/2, formula 12) for the 
low excitation regime 1] <C cjb- In regimes of stronger 
excitation, for higher electron-hole densities the screen- 
ing of interparticle Coulomb interactions in a semicon- 
ductor becomes important.^ Here we assume that the 
microscopic definition of J from Ref. 9 remains valid 
even in those regimes, but for calculation of J within 
the jellium model instead of the bare Coulomb inter- 
action potential and the respective Is-exciton hydro- 
genic wave function we use a proper screened potential 
and the corresponding ground-state eigenfunction of the 
Schrodinger equation. We utilize a natural, statically 
screened Coulomb interaction potential - the Yukawa po- 
tential:^ v{r) = e^(47r£o^B)~"'^ exp[— /^s^J/r, where r is 
the radius- vector modulus, and is the screening wave 
number. Accurate analytic approximations of the re- 
spective ground-state eigenvalue and eigenfunction are 
given, for example, in Ref. 20. Phenomena related with 
excitons take place at rather low temperatures, there- 
fore all the experiments on the SITm in semiconduc- 
tors were performed in this range of temperatures (e.g., 
8 K, like in Refs. 5 and 16]). That is why we use here 
the Thomas-Fermi form of the screening wave numbei^ 
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FIG. 4. (a) Numerical demonstration of a pulse breakup un- 
der propagation by using system ([5|. The electric field IV^I 
(solid), the polarization \p\ (dashed) and the inversion w (dot- 
dashed line) after a propagation of x = 60 are shown. The 
initial double pulse (dotted line) is '0(x = 0, T ) — 2'0(x = 
0, T) exp[i(^(x = 0,T)] in accordance with Eqs. (6b) and ([7|, 
while the parameters are A^ = 1 and J' — (similar to those 
at Fig. IsFa)). The areas Oj of electric field pulses are indi- 
cated, (b) Changes in the total IV^I -pulse area Oj with time 
after the propagation. The regions, which correspond to the 
areas of formed pulses, are indicated by the grey bars. 



= [/ie^(3A/x/7r)"^/^/(7r£o^B^^)]"^^^, where /i is the exci- 
ton reduced effective mass (it is 0.0969 tUq for CdSe),!^ 
and A/x is the total density of excitons (e-h pairs) gen- 
erated by an incident pulse. We estimate it here as 
Ax ~ I/{cnhuj^)^ I = cneo{hQ)'^ / {2dl^) is the incident 
pulse intensity.^ In this way we come to an approximate 
but natural and explicit relation between the model pa- 
rameter J and the incident field Rabi frequency - J can 
essentially decrease, when f]. Ax, and hv^ are large enough 
and the screening effects become significant. 

The limiting behavior of the SITm solitons is thereby 
straightforward. For the high excitation regime {Q ^ 
cjb) the screening is strong, electron- hole pairs are al- 
most decoupled, and J ^ 0, while Oj 27r in complete 
accordance with the results of analysis of microscopic 
SBEs limiting behavior in Ref. 11, For the low excita- 
tion regime (^7 <C cjb) the screening effects are negligible 
and J 13cc;b/3, while 6j becomes even smaller than 
unity, again in accordance with the results of analysis of 
the respective limiting behavior of SBEs.^ However, we 
should note that in the used model J remains finite and 
Oj ^ 1 even for Q below the low excitation regime, while 
in Ref. 10 in this case <C 1. This indicates a restric- 
tion of the used model - ^} should not be vanishing. Here 
and further in this section calculations are performed by 
using data on to = ^Jol = (r~^ + A^/4)~^/^ = 0.8 ps (r^ 
is the soliton temporal width) and A = 0, taken from the 
experiment on the SITm in CdSe.^ 

In the intermediate case ^ (for CdSe this cor- 
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responds to / ~ 750 MW/cm^), for which an approxi- 
mate doubhng of the incident field Rabi frequency should 
occur fISIIIl estimate J and Oj using typical incident 
pulse intensities employed in experiments on the SITm 
in CdSe:^^ / - 100 MW/cm^ According to the devel- 
oped formalism this yields o:^ 3.8 • lO^^cm"^, l//^s — 
8.94 nm, J 0.236 c^b, and hence Oj 1.087r (e.g., 
Oj e [0.82, 1.48]7r for / G [50,200] MW/cm^). Therefore, 
in good accordance with the known experiments^ for 
the considered intensities range the SITm effect indeed 
takes place for incident pulses of areas Oj tt. 

V. CONCLUSIONS 

In summary, a strong resonant MI effect in semicon- 
ductors has been demonstrated analytically on the basis 
of a known simple semiclassical macroscopic model ,1^^^ 



which originates from microscopic SBEs. New ways to 
shift the MI gain peaks, and therefore the generated fre- 
quencies have been discussed. The analytical results are 
completely confirmed by direct numerical simulations of 
the governing equations. The validity of the used macro- 
scopic model has been shown by studying the SITm 
effect in semiconductors, for which a good agreement 
with both existing theoretical^^ and experimentai^ 
results has been found. Further investigation of nonlin- 
ear light-semiconductor interplay beyond the SBEs for- 
malism (e.g., see Ref.[2T|) can lead to a revelation of novel 
aspects of the studied phenomena. 
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